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ABSTRACT 


This  thesis  investigates  the  effects  of  atmospheric 
turbulence  on  the  spatial  properties  of  the  phase  of  transmitted 
beams.   Experimental  measurements  in  the  form  of  a  series  of 
photographs  of  interference  patterns  between  two  beams  were 
made  using  a  modified  optical  heterodyne  system.   The  two  beams 
of  the  heterodyne  system  were  mixed  with  a  very  small,  control- 
lable angular  difference  between  their  directions  of  approach  to 
the  detector.   The  resulting  fringe  patterns  were  effected  by  the 
atmosphere  as  differences  between  succeeding  photographs  show. 
The  first  order  effect  of  the  atmosphere  was  found  to  be  a  tilt 
relative  to  the  direction  of  travel.   Tilt  differences  as  small  as 
5  seconds  of  arc  can  be  measured  using  this  system  with  assurance; 
smaller  values  could  be  confused  with  systematic  effects. 
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INTaODUCTION 

The  advent  of  the  laser  has  stimulated,  interest  in  communi- 
cation at  optical  frequencies.   The  potential  bandwidths  and  an- 
tenna gains  hold  great  promise,  but  the  degrading  effects  of  the 
atmosphere  on  propagating  optical  wavelength  beams  have  prevented 
immediate  widespread  usage.   Present  efforts  in  the  area  of  atmos- 
pheric optical  communication  are  mainly  directed  toxmrds  gaining 
an  understanding  of  the  atmospheric  phenomena  and  characterizing 
the  received  fields  after  propagation  through  the  atmosphere. 
This  information  is  necessary  to  build  a  statistical  model  of  the 
channel  at  optical  and  near  infrared  frequencies.   A  statistical 
model  will  allow  the  use  of  tools  from  modern  communication  theory 
to  establish  the  ultimate  limits  of  optical  communication  links 
utilizing  atmospheric  channels. 

Atmospheric  effects  on  propagated  fields  at  optical  and  near 

infrared  frequencies  can  be  broadly  put  into  t^ree  catagories! 

1.  Molecular  absorption 

2.  Particle  scattering 

3.  Turbulence 

The  first  two  are  more  or  less  suggestive  and  under  some  condi- 
tions these  two  effects  can  be  lumped  as  a  path-length  dependent 
loss  term.   This  assumption  depends  strongly  on  atmospheric  condi- 
tions and  propagating  wavelengths.   However,  with  a  clear  atmosphere 
this  model  is  reasonable  at  a  wavelength  of  6328/\.   This  thesis 
deals  with  the  third  effect,  turbulence. 

Turbulence  is  the  name  used  to  describe  the  effects  of  small 
variations  in  the  local  refractive  index  caused  by  turbulent  mix- 
ing of  the  atmosphere  under  wind  and  temperature  gradients.  At- 
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mospheric  turbulence  results  In  an  inhomogeneous  medium,  and  pro- 
pagation through  this  medium  causes  random  spatial  and  temporal 
fluctuations  in  amplitude  and  phase  of  the  received  field.   A 
"body  of  work  has  been  done  on  the  amplitude  fluctuations,  called 
scintillation,  with  beams.   There  has  been  less  emphasis  on  the 
phase  process  due  to  the  difficulty  of  phase  measurements. 

This  thesis  considers  a  method  of  phase  measurement  using 
an  optical  heterodyne  system  and.  evaluation  of  the  experimental 
measurements. 

Chapter  1  presents  a  brief  survey  of  the  pertinent  theory 
and  modeling  necessary  for  this  thesis. 

Chapter  2  describes  the  experimental  system  used  to  hetero- 
dyne at  6328.fi  and  the  data  analysis  system  in  general.   Appendix 
A  contains  a  more  detailed  explanation  of  the  system,  special 
devices,  and  electronic  circuits. 

"Chapter  3  presents  the  experimental  results.   Appendix  B  con- 
tains more  detailed  information  on  the  method  of  data  extraction. 

Chapter  4  considers  the  experimental  results  and  conclusions, 


CHAPTER! 

This  chapter  is  concerned  with  some  of  the  theoretical  re- 
sults developed  for  beam  propagation  through  the  turbulent  atmos- 
uhere  with  emphasis  on  those  results  Dertinent  to  the  experimental 
system  constructed. 
1-1  Atmospheric  Turbulence 

The  effects  of  atmospheric  turbulence  on  optical  wavelengths 
have  long  been  known  to  astronomers;  the  twinkling  of  starlight  on 
a  windy  night  is  a  common  example .   The  possible  use  of  the  laser 
for  atmospheric  communication  links  has  renewed,  interest  in  the 
study  of  atmospheric  turbulence  and  its  effects. 

Turbulence  is  small  scale  vatiations  of  the  index  of  re- 
fraction,  typically  on  the  order  of  1  part  in  10  ,  caused  by 
temperature  discontinuities  set  up  by  turbulent  mixing  of  the  at- 
mosphere.  These  refractive  index  changes  were  first  described  using 
a  Gaussian  correlation  function  which  was  an  analytically  convenient 
method  but  led  to  questionable  results.   A  second  path  taken  by 

Tatarski  applied  the  Kolmogorov  theory  of  turbulent  flow  at  very 

i 

large  Reynolds  numbers  to  the  index  variations.   This  analysis 

leads  to  a  second  moment  description  of  the  medium,  a  method  that 
is  not  alone  sufficient  to  describe  atmospheric  turbulence.   How- 
ever, the  fields  propagated,  through  the  lnhomogeneous  medium  will 
have  lognormal  statistics;  thus  a  second  moment  description  of  the 
medium  will  be  adequate  to  handle  the  received  field  statistics. 
The  results  of  Tatarski' s  work  show  that  the  turbulence  is  locally 
isotropic  and  homogeneous  and  has  a  structure  function  of  the  form 
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2  _      2/3 
Cn(R)p  Go#P£Lo 

Dn(p,R)=    <  (1) 

2  ,2/3,2 


Cn(R)^o        (fQ)  p<r<|o 


where 


Dn(p,R)=   [nC^)    -  nd^+p)]2  (2) 

p  ..   2r-.  -f  p       where  R  is  a  vector.  A   bar  will  be 
M   2  used,  to  denote  vectors. 

P  =  p 

2  _ 
C  (R)  =  constant  associated  with  the  magnitude  of  the 
n      turbulence  at  R 

n(r)  =  refractive  index  at  r 

*o  =  inner  scale  of  turbulence;  on  the  order  of  1mm 

and  related  to  the  breakup  of_laminar  flow.   At 

distances  greater  thanJ?o,  n(r)  tends  to  be  dis- 
continuous. 

Lo  =  outer  scale  of  turbulence;  on  the  order  of  meters 
and  related  to  the  gross  layer  structure  of  the 
atmosphere  or  the  distance  between  obstacles. 

The  model s  while  not  as  convenient  analytically,  is  generally 

description 
considered  better  than  the  Gaussian  correlation  function/ and  has 

been  the  basis  for  most  further  work. 

1-2  Propagation  through  the  turbulent  atmosphere 

The  available  theoretical  solutions  generally  are  for  infinite 

plane  wave,  sources  and  assume  a  "frozen"  atmosphere,  where  the 

turbulence  is  fixed  in  space.   The  frozen-in  structure  is  assumed 

to  move  across  the  transmission  path  due  to  the  component  of  wind 

normal  to  the  propagation  path  to  account  for  the  temporal  behavior. 

The  three  basic  approaches  to  the  turbulent  propagation  problem  are  I 

i 
using  a  Rytov  approximation  to  the  wave  equation,  solving  the  wave 
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2 

equation  for  the  mutual  coherence  function 5   and  a  geometric  optics 

3 
type  of  solution.    The  last  method,  geometrical  optics,  is  general- 
ly considered  inferior  since  it  does  not  account  for  diffraction 
effects  on  long  paths  Kith  small  scale  inhomogenetics. 
1-3  Received  Field  Parameters 


Using  the  "frozen"  atmosphere  assumption,  the  central  limit 
theorem,  and  one  of  the  solution  approaches,  it  can  be  shown  that 
a  time  sample  of  the  received  field  envelope  is  a  lognormal  random 
variable,  and  that  a  time  sample  of  the  received  field  phase  is  a 
Gaussian  random  variable.   There  is  evidence  to  expect  that  the 
form  of  these  statistics  is  unchanged  for  beam  propagation. 

Another  measure  of  importance  is  the  spatial  coherence  of  the 

phase  of  the  received  field.   Tatarski  has  derived  an  expression  for 

the  case  of  a  linearly  polarized  olane  wave  with  the  atmospheric 

2  _ 
turbulence  described  as  in  Eq.  1,  and  C  (R)  a  constant  independent 

_  1 

of  R.    Then 

225/3 
D<j>(p)  -  aCnk  Lp  (3) 

where  D^(p)  is  the  phase  structure  function,  k  the  wave  number  of 

the  radiation,  L  the  path  length,  and 

C   2.91  P  >  XL 

a  =    < 

Using  the  expression  for  DA(p")  in  Eq.  3,  a  coherence  diameter  for 


some  selected  mean-square  difference  in  phase  can  be  calculated. 

9  4 

One  example  quoted  in  Hovers ten  from  Goldstein   is 

"coherent  =  L?5  *  10~2  _xi  (4) 

CnL 
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Fried  has  used  a  form  of  Eq.  3  to  show  that  wavq^front  tilting  is 

5" 

the  largest  contributor  to  distortion. 

/mother  parameter  is  the  frequency  spectrum  of  the  amplitude 
and  phase  turbulence  processes.   Measurements  of  intensity  spec- 
trums  have  shovm  them  to  be  no  more  than  several  hundred  cycles 

V 

wide  for  most  conditions.   Measurements  for  the  phase  process  are 


not  as  plentiful,  but  recent  measurements  have  shown  that  the  spec- 

a 

trura  is  no  more  than  several  hundred  cycles  wide. 

X.r _il:  ^QjniUJlicAt iQJl __Model 

A   suggested  model  describing  the  effects  of  the  turbulent 
atmosphere  as  a  communication  channel, when  particle  scattering 
is  not  important,  is  a  multiplicative  channel  with  additive  back- 
ground  noise.   A  linearly  polarized  field  (chosen  for  simplicity) 
in  the  absence  of  losses,  turbulence  and  background  noise  would 
be  received  as 

Erec(r,t)  =  Real  JE(r,t)e    \  (5) 

When  turbulence  and  background  noise  are  present  the  received  field 
(again  in  the  absence  of  losses)  would  be 


E 
rec 


C   ~    J<»ttf(rtt)      _    J^i(r,t)  Ju)t? 
(r,t)  =  RealjE(r,t)e   e       +  E  (r,t)e        e   f 


(6) 


where 


e   *  '   represents  the  effect  of  the  turbulent  channel 

tf(r,t)  =  X(r,t)  +  J<jj(r,t) 
X(r,t)   represents  the  scintillation  process 
<£>(r,t)   represents  the  phase  distortion  process 
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E_^(T,t)e  represents  the  relevant  polarization 

component  of  the  background  radiation; 
a  complex  Gaussian  process  with  a  very 


wide  bandwidth  on  both  Enattd®n» 


^2^iJS§tl-_Ii§i§£2^ZGill£ 


Two  viaves  simultaneously  incident  on  a  surface  will  cause  an 
interference  pattern.   This  can  be  modeled  simply,  using  plane  waves 
of  the  same  frequency  and  linear  polarization.   When  the  two  waves 
beat  together  from  slightly  different  directions,  p_  and  fL,  a 
fringe  pattern  whose  period  is  related  to  the  difference  in  approach 
angle  is  created.   If  one  of  the  waves  has  been  sent  through  the 
atmosphere,  it  x\Till  have  a  phase  term  proportional  to  the  phase  dis- 
tortion process.   Thus  if  complex  notation  is  used,  suppressing 
the  e^  time  variations  for  simplicity,  the  fields  can  be  des- 
cribed as  (see  Fig.  1) 

-j0z  X(r,t)  34>(r,t)  (8) 

Erec(r,t)  =  E^e    e      e 

-JB^z-jP  x 


SLo(r,t)  =  E?e        *  (9) 


'Lo^'"'  ~  *2 


where 


-i 

h 


B?  I  =  //R2  .  «2      =s  —  =  ~I  and  g_  and  g0are 
Pz   *x  °   K       wave  vectors? 


Then  on  the  surface  (z=o)  shown  in  Fig.  1  the  total  field  can  be 
written  as 


X(p,t)  j<t>(p,t)      -jBvx 
ild       e        *  ^2 


=  S,e      e        +  E^e   x  (10) 


where 
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P  -  /xx  +  Zyy 


r  =   p  +  zzz 


Figure  1. 

Optical   Heterodyne 

Geometry 


Then  the  intensity  on  the  surface  is 


I(p,t) 


1     f»  F*  1 

^    |^total~totalj 


JL 

2rn 


2    2X(p,t) 
e  + 


hJ 


+    2 


Wl  M 


X(p,t) 


e 


cos[Bxx-cp(p,t)-a]J 


(11) 


where -Y)    is  the  impedance  of  the  atmosphere  and  a  represents  the 
phase  angle  between  the  complex  amplitudes  EL  and  E?.   Thus  a  por- 
tion of  the  intensity  is  a  term  phase  modulated  by  the  atmospheric 
phase  distortion.   In  the  experimental  section  of  this  thesis  a 
system  to  exploit  this  effect  is  described  in  detail. 
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CHAPTER  2 

This  chapter  describes  the  system  developed  to  create  a 
spatial  enterference  pattern  between  two  beams  of  coherent  light, 
one  modulated  by  atmospheric  turbulence. 
2  -  1  System  Layout  and  Eq uipment 

The  basic  arrangement  for  heterodyning  is  shown  in  Fig.  2. 
The  source  was  a  Spectra-Physics  Model  119  Helium-Neon  gas  laser- 
operating  at  a  wavelength  of  6328A.   The  Model  119  laser  has  a 
rated  continuous  output  power  of  0.1  milliwatt  although  measure- 
ments indicated  slightly  higher  values  during  experimentation.   The 
beam  from  one  end  was  collimated  using  a  Spectra-Physics  Model  331- 
333  spatial  f ilter-collimator  and  transmitted  over  an  atmospheric 
path.   The  down  range  reflector,  shown  in  Fig.  3  and  Fig.  4,  is  a 
sixteen-inch  diameter  flat  mirror  mounted  in  an  aluminum  collar  and 
frame.   The  collar  and  frame  have  precision  adjustments  through  screw- 
arm  systems  that  allow  very  accurate  pointing  of  the  return  beam. 
The  beam  from  the  opposite  end  of  the  laser  was  reflected  off  two 
flat  first  surface  mirrors  mounted  on  precision  adjustment  holders 
to  reach  the  combining  beam  splitter.   No  lens  was  used  in  this  non- 
atmospheric  path  to  minimize  extraneous  lens  aberrations.   As  this 
beam  was  close  to  plane  wave  conditions  by  the  time  it  reached  the 
combining  beam  splitter,  its  curvature  was  easily  matched  by  the 
telescope  output  using  the  telescope  eyepiece. 

The  receiving  telescope  is  a  three-inch  reflector  with  an 
eighteen-inch  focal  length.   The  lens  from  a  Praktica  FX-3  camera 
(f/2.8,  5°  millimeter  focal  length)  was  used  as  the  eyepiece  for 
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several  reasons;  the  50  millimeter  focal  length  meant  the  telescope 
output  beam  size  closely  matched  the  local  oscillator  beam  size 
and  the  focusing  ring  of  the  eyepiece  allowed  adjustments  to  be 
made  with  a  minimum  of  mechanical  movement. 

The  two  signals  were  combined  using  a  fifteen-millimeter 
prismatic  beam  splitter  mounted  on  a  fine  adjustment  holder.   This 
mount  was  used  to  manipulate  the  fringe  spacing  and  orientation 
since  relatively  small  changes  in  angular  separation  can  be  made 
most  easily  at  this  point. 

All  of  this  equipment,  with  the  exception  of  the  large  doxm- 
range  reflector9  was  bolted  down  on  a  two-inch  thick  aluminum  op- 
tical bench  housed  in  a  penthouse  on  the  roof  of  Building  20  at  H»I.T. 
(See  Figs  5  and  6).   The  downrange  reflector  was  housed  in  a 
protective  box  on  the  other  end  of  the  roof  of  Building  20,  mak- 
ing the  atmospheric  path  approximately  2  00  meters  round  trip. 
2-2  Pattern  Detection 


The  interference  pattern  of  the  two  beams  was  recorded  on 
black  and  white  film  using  a  Praktica  FX~3  35^m  single  lens  reflex 
camera.   The  351TU^  negatives  were  developed  using  Kodak  formula  D-?6 
developer,  standard  for  the  Tri-X  film  used,  and  were  mounted  in 
holders  to  be  used  as  slides  for  later  projection.   The  negatives 
were  used  directly  for  this  projection  step  since  the  pattern  of 
light  and  dark  fringes  was  the  same  whether  or  not  the  image  inten- 
sities had  been  reversed. 

The  fringe  patterns  were  examined  in  detail  using  a  slide 
projector  to  put  the  image  on  a  screen.   This  image  was  then 
scanned  using  a  Packard  Bell  Model  920  television  camera.   The 
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video  output  of  the  television  camera  was  displayed  on  a  Tech* 

tronix  Model  $k>$   oscilloscope  as  described  in  section  2-*K 

^.J^XJ5^§Zi£lP2l5ilL.PX9P^eil^£ 


The  alignment  between  the  source  and  receiver  was  accomplished 
by  finding  the  laser  (not  turned  on)  image  in  the  receiving  tele- 
scope and  down-range  mirror  system  during  daylight  hours.   Thus 
the  optical  axis  between  laser,  mirror,  and  telescope  was  estab- 
lished.  The  laser  beam  could  by  angular  adjustment  of  the  laser 
mount  then  be  put  on  the  mirror  at  night  to  have  this  beam  direction 
coincide  with  the  optical  axis. 

At  this  point  the  beam  emerging  from  the  opposite  end  of  the 
laser  was  positioned  with  the  small  flat  mirrors  so  that  it  pointed 
at  the  center  of  the  beam  splitter.   All  further  adjustments  were 
made  using  the  beam  splitter.   The  local  oscillator  signal,  making 
a  right  angle  change  in  direction  was  much  more  sensitive  to  the 
movement  of  the  beam  splitter  than  the  beam  coming  out  of  the  tele- 
scope  going  straight  through.   Thus  the  local  oscillator  beam  could 
be  moved  without  disturbing  the  telescope  output. 

The  beams  coming  out  of  the  beam  splitter  could  be  safely  ob- 
served by  eye  after  they  had  passed  through  a  Kodak  ■$  2   Wratten 
filter.   Without  an  eyepiece  the  beams  appeared  as  pencil  points 
which  could  be  adjusted  to  coincide.   An  eyepiece  could  then  be 
used  to  see  the  resulting  fringe  pattern. 

After  establishing  a  fringe  pattern,  the  filter  was  removed 
and  the  Praktica  FX-3  camera  without  lens  was  used  to  record  se- 
quences of  three  to  six  consecutive  frames  of  the  same  fringe  pat- 
tern ie,,  three  to  six  frames  were  exposed  without  altering  the 
optical  system. 
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2  -  ^  Data  Pre sen tat i on 

The  slides  used  for  analysis  were  scanned  using  a  Packard 
Bell  920  television  camera  synchronized  to  an  external  clock. 
This  same  clock  output  was  used  to  trigger  the  oscilloscope. 

The  video  output  was  displayed  in  two  distinctive  manners. 
The  first  method  was  designed  to  obtain  an  overall  picture  of  the 
fringe  process  across  the  aperture,  the  second  method  to  look  at 
narrow  strips  of  the  fringe  process  across  the  aperture. 

The  first  method  is  shown  in  block  diagram  form  in  Figure  7. 
■The  horizontal  sweep  lines,  modulated  by  the  alternating  dark  and 
light  strips  of  the  projected  fringe  patterns,  were  displayed  one 
on  top  of  the  other  by  setting  the  time  base  of  the  oscilloscope 
to  the  horizontal  sweep  time  of  the  television  camera.   The  ver- 
tical deflections  of  each  line  represent  the  tones  of  grey  as 
scanned.   Then  a  ramp  signal  was  added  to  the  horizontal  sweep 
lines  to  separate  them  out  vertically.   This  ramp  signal  was 
synched  to  the  vertical  sweep  time  of  the  television  camera. 
Thus  a  quasi -television  picture  is  made  with  vertical  voltage  de- 
flections replacing  shades  of  grey. 

The  second  method  used  to  display  the  television  scan  was  to 
select  out  one  individual  horizontal  scan  line  with  the  system  as 
shown  in  Figure  8.   The  "/Jclock  delayed  by  B  clock"  mode  of  the 
Techtronix  5^5  oscilloscope  was  used.   The  A   delayed  by  B  mode  is 
a  system  where  first  a  ramp  signal  set*  to  the  B  clock  time  base 
is  generated.   This  ramp  signal  is  compared  to  a  value  dialed  in 
by  the  scope  operator.   When  the  two  signals  are  equal,  the  A  clock 
is  armed,  to  take  the  next  trigger  it  sees  and, when  triggered, to 
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Figure  7.   Full  Aperture  Display 
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sweep  once  on  its  time  base.   Thus  one  sweep  (A)  In  a  sequence  (3) 
can  be  selected.   By  having  the  B  cloclc  synched  to  the  vertical 
sweep  of  the  television  camera  and  displaying  with  the  A  clock   . . 
synched  to  the  horizontal  sweep  of  the  television  camera,  a  select- 
ed horizontal  sweep  of  the  television  camera  can  be  displayed. 
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CHAPTER  J3. 

Chapter  3  describes  the  experimentally  observed  results  and 
the  descriptions  of  the  phase  processes  inferred  from  this  data. 

The  individual  intensities  of  the  local  oscillator  beam  and 
the  transmitted  beam  after  propagation  through  the  atmosphere  and 
receiver  optics  are  shown  in  Figures  10  and  11  respectively. 
These  photographs  were  made  by  simply  blanking  off  the  other  beam, 
and  so  represent  typical  intensity  distributions  at  the  film  sur- 
face for  each  of  them.   The  picture  of  the  intensity  over  the 
receiver  aperture  shows  the  effects  of  scintillation,  ie.,  a  non- 
uniform distribution  of  intensity  with  a  spatially  random  pattern. 
The  blank  spot  in  the  middle  and  the  strip  running  off  to  the 
right  in  Figure  11  are  due  to  the  secondary  mirror  of  the  Newton- 
ian telescope  used  as  the  receiver.   Although  this  mirror  blanks 
off  an  appreciable  portion  of  the  photograph,  enough  of  the  aper- 
ture is  seen  to  allow  analysis. 

When  the  two  beams  are  present  simultaneously  and  allowed 
to  interfere  on  the  film  surface,  the  fringe  patterns  shown  in 
Figures  12,  13  and  14  are  created.   It  was  found  through  experi- 
mentation that  easily  discernable  changes  in  the  fringe  spacing 
over  the  aperture  could  not  be  seen  unless  the  number  of  fringes 
was  reduced  to  approximately  twenty  or  less.   If  the  model  of  Ap- 
pendix A  is  assumed,  then  equation  (A-12)  describes  the  fringe 
pattern  as 
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Figure  10.   Local  Oscillator  Beam  Intensity  at  Detector 
Surface 


Figure  11.   Transmitted  Beam  Intensity  at  Detector  Surface 
after  traversing  Atmospheric  Path  and  Receiving 
Optics. 
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i    C     ,   ?-  2  2   2,         2    2X(p") 

«P-)    -  ^fM       +  K     +   u   L   jHj     e 

i       ,  ,  X<P>  -     -       J  ") 

+   2WLJUJ    ]   U2je  cos[Bp'  p   -<}>(  p)-a]  V  (A-12) 

where 

t^  represents  the  amplitude  of  the  envelope  of  the  field 
sent  through  the  atmosphere 

U2  represents  the  amplitude  of  -the  envelope  of  the  local 
oscillator  field  where  a  is  the  phase  angle  between  E, 
2   and  Ep  ~"1 

N  represents  the  average  collected  background  radiation 
intensity 

L  represents  the  losses  in  sending  a  beam  through  the 
atmospheric  channel  and  in  collection 

B  represents  the  direction  cosine  between  the  two  waves 
at  the  detector  surface  (see  Fig.  1) 

rj   represents  the  impedence  of  the  atmosphere 

M  represents  the  magnification  of  the  telescope  and  eye- 
piece 

represents  the  scintillation  process  due  to  turbulence 
(£>(p)  represents  the  phase  distortion  process  due  to  turbulence 
This  intensity  expression  reflects  the  fact  that  this  pattern  is 
assumed  to  be  a  time  sample  of  the  fields.   It  contains  only  spa- 
tially varying  terms  and  implicitly  assumes  that  this  pattern  is 
taken  for  a  particular  time  t.   At  some  later  time  a  photograph 
would  have  another  time  evoltion  of  the  sample  functions  X(p) 
and  <p( p). 

As  shown  in  Appendix  B,  the  geometry  can  be  selected  so  that 
the  fringes  occur  along  one  of  the  axes  in  the  p"  plane.   With  this 
geometry  the  direction  cosine s  B  ,  can  be  related  to  the  difference 
in  approach  angles, 2  9    of  the  two  beams  so  that  the  intensity  be- 
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Figure  12.  Sequence  A  Fringe  Patterns.  These  consecut 
photographs  will  be  analyzed  in  the  followin 
figures:   Figs.  15 

18 

22 

25 
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Figure  13.   Sequence  B  Fringe 
Patterns.   These  three  con- 
secutive photographs  will  be 
analyzed  in  the  following 
figures! 


Figs. 


16 
19 
23 
25 
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Figure  14. 


Sequence  C  Fringe  Patterns.   These  consecutive 
photographs  will  be  analyzed  in  the  following 
figures:   Figs.  17 

20 

24 

27 


comes 
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2  2  2    2  2    2X(p) 

+  N     +  u  L  1  UJ     e 


X(P)  ^ 


+  2uL(%||H2|e  oosC^S  -(|>(p)-a]t  (12) 


where 


a  =  the  phase  angle  between  E-,  and  E? 
sin£  =   |Bp|  /J6|  £r£ 
and  £4C1 

It  can  be  seen  that  the  fringes  are  caused  by  the  phase  modulated 
sinusoid  riding  on  the  intensities  of  the  individual  signals. 
Further,  it  can  be  seen  that  if  the  period  of  the  sinusoid,  £■, 
is  small  compared  to  the  distance  over  which  significant  changes 
occur  in  the  values  ofa>(p),  the  phase  distortion,  individual 
fringe  lines  will  not  experience  much  perturbation.   Since  a 
maximum  of  twenty  lines  distributed  over  the  three  inches  of 
aperture  show  appreciable  perturbations  of  the  fringe  lines,  a 
minimum  significant  distance  ln<p(*p)  can  be  estimated  at  (3/20) 
inches  =  k   millimeters.   This  rough  estimate  of  significant  di- 
mensions in(p(p)  conforms  to  data  given  in  the  literature  and 
general  experience  from  longer  paths. 

The  amplitude  scintillation  in  the  three  sets  of  fringe  pat- 
terns does  not  appear  to  be  severe,  but  this  is  due  to  a  selection 
process  for  just  this  quality. 

Another  notable  feature  of  these  slides  is  the  spatially  vary- 
ing fringe  pattern  across  the  aperture.   The  fringe  spacing  in 
Figs.  12  and  13  are  narrower  on  the  left  side  of  the  aperture. 
Especially  in  the  photographs  of  Fig.  13,  one  can  see  that  the 
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j  . 

fringe  spacing  not  only  varies  spatially,  but  also  chances  radi- 
cally from  one  photograph  to  the  next.   The  photographs  in  Fig.  Ik 
show  a  rotation  of  lines  between  photographs.   The  photographs  of 
Fig.  12  were  taken  just  prior  to  those  in  Fig.  13,  while  those  in 
Fig.  1^  were  taken  on  another  night  without  any  appreciable  modi- 
fication of  the  system  having  taken  place. 

A  photograph  of  the  system  through  a  non-atmosDheric  path 
could  not  be  obtained  unless  the  system  as  installed  was  changed 
significantly.   The  divergence  of  the  wave  emerging  from  the  laser 
or  laser— collimator  when  beamed  into  the  receiver  telescope  was 
such  that  its  focusing  point  was  moved  too  far  out  to  use  the  ca- 
mera eyepiece  as  mounted.   This  camera  eyepiece  was  needed  to  match 
curvatures  with  the  LOwave.   Other  eyepieces  required  optics  in 
thelLO  path.   When  the  divergence  was  reduced  so  that  the  system 
could  be  used  "as  is%  the  aperture  was  not  fully  illuminated,  and 
very  little  of  the  aperture  could  be  seen  around  the  sedondary  mir- 
ror.  Thus  system  effects  can  not  be  directly  taken  out.   However, 
the  availability  of  consecutive  photographs  allows  analysis  of  the 
changes  in  successive  photographs.   The  ootical  system  is  the  same 
between  photographs  except  for  advancing  the  film  in  the  camera, 
so  that  changes  in  the  fringe  spacing  would  indicate  atmospheric 
effects.   For  this  reason,  the  photographs  of  Fig.  1^4-  were  included 
to  obtain  another  consecutive  look  $.t  the  lumped  system-atmospheric 
effects  on  a  night  different  from  that  of  Figs.  12  and   13. 
3,.,-  1  Full  Aoerture  Analysis 

In  an  effort  to  obtain  a  feel  for  the  processes  that  occur 
over  the  entire  aperture,  the  slides  were  projected  on  a  screen 
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and  scanned  with  a  television  camera.   The  output  of  the  television 

camera  was  displayed  on  an  oscilloscope  using  the  first  of  the 
methods  described  in  section  2-4.   This  method  is  essentially 
the  construction  of  a  television  picture  but  with  the  shades  of 
grey  replaced  by  vertical  deflections  on  the  individual  scan  lines. 
This  method  gives  the  appearance  of  a  three  dimensional  map  of  the 
intensity  over  the  aperture  and  the  results  are  shown  in  Figs.  15, 
16  and  17. 

The  main  interest  in  displaying  in  this  manner  is  to  get 
some  feel  for 

(1)  Phase  changes  along  the  fringes 

(2)  Phase  changes  across  the  fringes 

(3)  Scintillation 

(4)  Systematic  effects 

The  deep  troughs  on  the  edges  of  the  photographs  in  Figs. 
15,  16  and  1?  are  the  locations  of  black  non-reflecting  tape  put 
on  the  projection  screen  to'  aid  in  centering  the  picture.   In  e- 
yery  photograph,  the  fringe  spacing  increases  towards  the  right. 
This  is  considered  a  system  constant  and  as  such  will  be  dis- 
counted in  further  analysis. 

The  fringe  lines  tend  to  curl  towards  the  lower  right-hand 
corner.   This  curl  is  stronger  in  the  lower  half  of  the  aperture 
and  can  be  seen  particularly  well  in  photograph  of  Fig.  17.   This 
Is  an  example  of  the  phase  distortion  along  the  f^ringe  lines,  A 
slightly  different  effect  is  the  amplitude  modulation  along  the 
fringe  lines  caused  by  phase  distortion. 

The  scintillation  can  be  seen  best  on  the  peak  lines  as  on 
amplitude  modulation.   Where  the  modulation  produced  by  the  fringes 
is  shallow,  generally  on  the  right  hand  sides  of  the  photographs, 
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Figure  15. 


Sequence  A  Full  Aperture  Display  of  the 
Fringe  Pattern  of  Photograph  #1  andd*  2 
by  Method  of  Section  3-2. 
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Figure  16. 


Sequence  B  Full  Aperture  Display  of  the 
Fringe  Pattern  of  Photography  1  by  Method 
of  Section  3-2. 


Figure  17. 


Sequence  C  Full  Aperture  Display  of  the 
Fringe  Pattern  of  Photograph  # 1  by  Method 
of  Section  3-2 
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5cintillation  tends  to  swamp  out  the  detail  of  t!he  pattern. 

The  analysis  of  the  phase  process  across  the  fringe  grain 
illl  be  done  in  detail  in  the  next  two  sections.   It  is  worth- 
while noting  at  this  ooint  that  the  fringe  pattern  is  generally 
/isible  and  enough  definition  exists  to  justify  a  look  in  detail. 

J.Trr  3   Phase  Analysis 

Using  equation  (B-2)  as  the  model  for  this  fringe  line  sys- 
tem, the  intensity  is 

-,  r  2     2     2  2     2  2X(p") 

I(p>  =^\\M   +  N  +M  L  N  e 


+   2mL|u2|  I  Uj  e  cox[^^  -  <|>(x9yo)~a]j      (12) 

If  a  narrow  strip  in  the  y  dimension  is  examined,  narrow  enough 
so  that  the  processes  do  not  change  with  y,  the  intensity  ex- 
pression becomes 

,  f    2     2    2  2     2  2X(x,yn) 
Kx,y0)  =  Yi\   (|^2|   *  N  *  u  L  |%|  e 

I'  +  2nLJy2j  JuJ  e   '  °  cos[^^  -  (j>(x9yo)-ai      (13) 

Assuming  that  the  spatial  sinusoid  is  visible  despite  the  random 
amplitude  modulation  from  scintillation,  the  distance  between  suc- 
cessive peaks  or  valleys  can  be  measured.   This  indicates  the  sin- 
usoid argument  going  through  2tt  radians  or  one  period.   Thus, 
choosing  the  peak  lines  as  the  point  to  measure  period  length  on, 

we  know  over  one  period 

l2tt£x  _^fv   „  N1  _  r2rr£x 

where 


[£!!*&  -.(h(x,yo)]  -  [^  -Cp(x»,yJ]  =  2tt  (1*0 


x  and  x'  are  the  points  where  two  adjacent  maximum  values 
of  the  sinusoid  are  located.   In  the  remainder 
of  this  thesis  x  and  x'  will  Indicate  adjacent 
maxima. 
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Sincc   x  and  x'    can  be  measured,    we   get 

<j>U,y0)  -<jx*\y0)  =  2tt[|(x-x»)  -  i]  (15) 

Ihus  over  the  distance  (x-x»)  we  get  a  measure  of  how  the  phase 
Listortion  process  acted.   If  we  continue  this  analytic  procedure 
across  the  fringe  pattern,  vie  have  some  notion  of  how  <b(x,y  ) 
changes  with  x. 

It  has  been  shown  by  Fried  that  geometrical  interpretations 
fen  be  placed  on  the  form  of  (J)(x,yo).   If  a  Taylor  type  expan- 
sion is  made  of  &)(x9y    ),  we  get 

T  O 


$ 


(x,yQ)   =  A  +  Bx  +  Cx2  +   .    .    .  (16) 


+   2uL 


Cutting   this   expression  for    Cp(x,y    )    into  Eq.    13  gives 

if;      ,2  2  2    2  2    2X(xsy    ) 

Kx,y0)    =  ^JJuJ       +  N     +  u  L   l^j    e 

.  ,      .    X(x,y   )  9    n  9 

-1    |~2   e  cosD^pc  -   a  -   A   -   3x   -   Cx    .    . 

(17) 

Tere  it  can  be  seen  readily  that 

A  represents  a  phase  origin  relative  to  the  position 
coordinates  x 

Bx  represents  a  linear  phase  change,  a  tilt  relative  to 
the  direction  of  travel 

2 
Cx   ...  represents  a  spherical  phase  change,  representing 

a  spherical  wavefront  on  the  received  wave 
tnd  the  coefficients  give  a  good  measure  of  the  importance  of  each 
term;  the  number  of  terms  needed  to  adequately  represent  the  phase 
>rocess  indicate  in  some  sense  how  simple  the  phase  changes  are. 
There  are  two  problems  involved  with  applying  this  analysis 
.irectly  to  the  fringe  patterns  made  for  this  thesis.   The  first 
s  that  no  non-atmospheric  pattern  exists  to  take  out  the  purely 
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systeraatic  effects.   The  only  method  of  removing  these  systematic 

effects  is  to  examine  consecutive  patterns  and  remove  the  constant 
affects;  this  leads  to  the  second  problem.   The  only  distinctive 
feature  that  can  be  used  for  locating  the  same  narrow  strip  in 
consecutive  patterns  is  the  secondary  mirror.   Thus  without  taking 
(Lata  for  several  entire  slides,  which  would  take  a  prohibitive 
mount  of  time,  only  limited  strips  can  be  handled  in  each  pattern 
or  comparison  with  the  next.   Fortunately,  the  regions  just  above 
md  below  the  secondary  mirror  are  not  unique  compared  to  the 
emainder  of  the  fringe  patterns  in  the  rest  of  the  aperture,  as 
jan  be  seen  by  examining  Figs.  12,  13  and  1*+.   Thus  these  regions 
ire  a  reasonable  choice  for  analysis. 
|-  k  Narrow  Strip  Anal y s is. 

To  examine  how  the  fringe  patterns  made  for  this  thesis 
[hanged  in  phase,  narrow  strips  above  and  below  the  secondary  mir- 
or  were  measured  in  detail.   The  slides  containing  the  fringe  pat- 
ems  were  urojectetf  on  a  screen,  then  scanned  using  the  televi- 
:ion  camera.   The  fringes  were  aligned  so  that  the  horizontal 
can  of  the  television  camera  was  across  the  grain  of  the  fringe 
>attern.   Selected  horizontal  lines  from  the  television  camera 
rere  displayed  and  the  locations  of  the  peaks  of  the  fringe  lines 
[ere  measured  for  five  to  fifteen  successive  horizontal  lines. 
Ms  corresponds  to  a  strip  one  to  three  millimeters  wide  over 
;he  objective  aperture.   Typical  displays  of  these  horizontal  lines 
Te   shown  in  Figs.  18,  19  and  20,  where  the  top  line  corresponds 
o  a  scan  line  in  the  top  half  of  the  fringe  pattern;  the  middle 
.ine  to  a  scan  line  In  the  middle,  and  the  bottom  line  to  a  scan 
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Figure  18.   Sequence  A 
Typical  Horizontal  Lines 
from  Television  Camera 
for  Photograph  ^ 1 


Figure  19.   Sequence  B 
Typical  Horizontal  Lines 
from  Television  Camera 
for  Photograph  **  1 
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Figure  20.   Sequence  C 
Typical  Horizontal  Lines 
from  Television  Camera 
for  Photograph  $  1 


VV-w 
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line in  the  bottom  half  of  the  fringe  pattern.   The  relatively 
small  modulation  in  the  middle  of  the  middle  lines  is  due  to  the 
secondary  mirror.   The  relatively  deep  excursion  of  the  scan  lines 
at  the  edges  is  caused  by  black  non-reflecting  tape  placed  on  the 
screen  to  provide  a  means  of  centering  the  fringe  pattern  on  the 
television  monitor.   The  spatial  scale  of  the  oscilloscope  display 
can  be  related  directly  to  the  spatial  scale  of  the  aperture.   The 
conversion  is  0.1*0  inches  of  aperture  is  equal  to  1  centimeter  of 
oscilloscope  display.   This  conversion  will  be  applied  to  all  fu- 
ture figures  that  deal  with  the  aperture. 

Typical  unaveraged  strips  are  shown  in  Fig.  21,  where  the 
peak  lines  are  drawn  in.   The  locations  of  these  lines  are  simply 
averaged  along  the  fringes.   The  distance  between  the  averaged  lo- 
cations of  the  peak  lines  across  the  aperture  is  plotted  in  Figs. 
22,  23  and  24.   From  Eq.  15  this  distance  between  peak  lines  can  be 
related  to  the  phase  distortion  and  fringe  period  as 

j    (i-x«)  =  F^~  [<j?(x,yo)  -<t>(x»,y0)  +  2tt]  (18) 

Thus  measuring  (x-x')  gives  a  number  proportional  to  <p  over 
the  distance  (x-x*).   Further  geometric  interpretations  about  the 
form  of  &  can  be  made  as  successive  (x-x')  values  are  measured  and 
plotted  across  the  aperture.   A  direct  measure  of  the  spatial 
frequency  ™  is  not  available,  however  this  term  can  be  assumed 

A 

constant  across  the  aperture.   Then  any  changes  in  the  (x-x') 
across  the  aperture  can  be  attributed  to  changes  in  <P  (x,yo)  and 
any  changes  in  (x~x')  between  photographs  can  also  be  attributed 
to  the  change  in  (j)(x,y  )  with  time. 

This  method  of  plotting  the  average  distance  between  peaks 
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;  Figure  21.   Typical  Unave raged  Locations  of  the  Peaks  of  a 

Fringe  Pattern.   This  example  comes  from  Photo- 
graph 1  of  Sequence  B. 
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Figure  22.   Average  Spatial  Period,  of  Fringe  Pattern 
over  Objective  Aperture   Sequence  A 
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Figure   23.      Average   Spatial   Period    of   Fringe   Pattern 
over   Objective   Aperture      Sequence   3 
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Figure  2k,      Average  Spatial  Period  of  Fringe  Pattern 
over  Objective  Aperture   Sequence  C 


gives  a  visual  display  of  the  general  type  curvature  of  the  phase 

process  across  the  aperture  and  seems  to  show  a  general  quadratic 

form  indicative  of  spherical  distortion  of  the  incoming  wave. 

Another  method  of  plotting  is  to  use  the  averaged  location  of 

the  peak  lines  on  the  abcissa  and  to  plot  the  2rr  radian  change  of 

the  sinusoid  argument  between  successive  peaks  on  the  ordinate. 

This  is  shovm  in  Figs.  25?  26  and  27.   These  graphs  have  a  linear 

mean  square  estimate  line  plotted  for  reference.   The  first  aspect 

to  note  is  the  dominance  of  the  linear  term  in  every  plot.   The 

lines  have  some  higher  order  deformations,  "but  the  linear  term  is 

obviously  most  important.   The  second  aspect  to  note  is  the  change 

in  slope  between  photographs  in  Fig.  26.   This  means  a  significant 

change  in  tilt  of  the  incoming  wave  has  taken  place.   We  can 

equate  the  argument  of  the  sinusoid  in  Eq.  13  to  the  estimate  lines 

in  Figs.  25,  26  and  27  as 
2ttjL  x 
r-y—  -<£(x,y0)-a]  -  kxx  +  kg  (19) 

where 

v  _  ™2JL   the  linear  estimate  line  slope  measured  for 
1  '    2  each  photograph 
x 

|j,,    the  covariance  between  the  averaged  location  of  the 
peak  lines,  x,  and  the  total  phase  change,  y,  at 
that  x 
2 

o     the  variance  of  the  averaged  location  of  the  peak  lines 

m   the  mean  of  the  averaged  locations  of  the  peak  lines 

m   the  mean  of  the  total  ohase  change  across  the  auerture 

ky       ull 
kp  =  m  -  ——-     m  the  linear  estimate  intercept  for  a 
y   c       photograph 

The  change  of  this  argument  between  successive  photographs  is 

^(x,yQ)  -  (j}'(x',yQ)   =  (l^-k'^x  -  (k2-k«2)  (20) 
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Figure  2$.      Total  Phase  Change  Across  the  Receiving 
Aperture  Secuence  A 
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Figure  26.   Total  Phase  Change  Across  the  Receiving 
Aperture   Sequence  B 
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'igure  26.   Total  Phase  Change  Across  the  Receiving 
Aperture   Sequence  B 
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Figure  27e   Total  Phase  Change  Across  the  Receiving 
Aperture  Sequence  C 
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where  the  primes  refer  to  the  succeeding  photograph.   The  (k?-k'p) 
is  not  meaningful  since  many  steps  in  the  analysis  introduce  arbi- 
trary references  in  phase.   However,  the  change  in  slope  is  easily 
measured.   These  changes  in  k,  are  shown  in  Table  1  and  the  equiva- 
lent change  in  tilt  of  the  incoming  wave  is  also  listed. 

^k(  radians/inch) 

Above  secondary  mirror      1.2rr 
Figure  25 

Below  secondary  mirror      2n 

Above  secondary  mirror  k.  25",T/.  75^ 
Figure  26 

Below  secondary  mirror  8  e  7517 

Above  secondary  mirror        .15tt 
Figure  2? 

Below  secondary  mirror        .  75TT 

Table  1    Measured  Changes  in  Tilt  in  Three  Gets  of  Consecutive 
Photographs  (from  Figs.  12,  13  and  1*0 

As  is  discussed  in  Appendix  A,  tilt  changes  less  than  15 
microradians  (3*2  seconds  of  arc)  will  be  indistinguishable  from 
possible  mirror  tilt  caused  by  vibrations.   This  limit  is  exceeded 
jonly  in  Fig.  26.   If  the  photographs  used  to  generate  Fig.  26  are 
[reexamined  (see  Fig.  13) s  it  can  be  clearly  seen  that  the  fringe 
jspacing  changes  dramatically  from  the  first  to  second  photograph, 
an  effect  not  seen  in  Figs.  12  and  1*4-. 

In  an  effort  to  determine  the  higher  order  deformations,  a 
log-log  plot  was  constructed  using  the  average  location  of  peak 
lines  as  the  abscissa  and  the  phase  changes, with  the  linear  esti- 
mate subtracted  out,  as  the  ordinate.   However,  only  systematic 
effects  are  visible. 
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This  chapter  summarizes  the  results  presented  in  Chaper  3? 
and  draws  some  conclusions  about  the  general  procedure.   Some  sug- 
gestions for  improvement  of  this  system  are  also  made. 

Using  the  methods  outlined  in  Chapter  3,  it  is  clearly  pos- 
sible to  measure  phase  in  a  narrow  strip  across  the  objective 
aperture.   Using  the  resulting  phase  plots,  it  is  possible  to  com- 
pute a  change  in  tilt  between  successive  photographs.   In  this  ex- 
periment, the  linear  terms  and  systematic  effects  masked  possible 
higher  order  deformations  caused  by  atmospheric  turbulence,  but 
with  more  care  on  a  longer  path  these  effects  should  be  measurable. 

The  details  of  the  analysis  show  that  the  minimum  tilt  changes, 
1  separate  from  mirror  induced  tilt,  measurable  by  the  present  sys- 
tem was  approximately  5  seconds  of  arc.   It  is  interesting  to  note 
that  the  measured  tilt  ranged  as  high  as  22,5  seconds  of  arc,  a 
large  value  for  a  three  inch  aperture  and  a  short  path. 

When  the  fringe  pattern  produced  by  the  two  beams  is  looked 
at  through  a  Wratten  filter,  the  system  dynamics  can  be  seen  graph- 
ically.  The  fringe  lines  waver  along  their  length  and,  a^s  their 
spacing  gets  wider,  a  bowing  travelling  across  the  pattern  is  of- 
ten seen,   A  third  effect  is  a  shimmering  of  the  entire  pattern, 
possibly  indicating  a  time  change  of  phase  origin  suggested  in 
section  3-3* 

The  system  of  oscilloscope  display  and  visual  reading  is 
clearly  a  tedious  method  of  data  taking.   The  amount  of  data 
that  can  be  taken  in  this  way  is  small,  and  is  not  compatible 
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with  machine  computation.   Judgement  can  be  exercised  to  pick 
out  the  peak  lines  from  scintillation  when  this  display  method  is 
used ,  but  the  peaks  and  valleys  of  the  fringe  pattern  are  the  only 
parts  that  can  be  used  for  visual  data  taking. 

The  lack  of  a  system-only  series  of  photographs  was  a  hin- 
drance to  the  analysis.   Some  further  feel  for  the  motions  of  the 
down- range  reflector  would  have  been  invaluable  in  interpretating 
the  results. 
4  _-.,.2_ Suggestions  for  Improvement 

The  first  suggestion  would  be  to  set  up  a  path  in  the  la- 
boratory long  enough  so  that  the  divergence  of  the  wave  entering 
the  telescope  is  not  too  fast.   At  the  same  time,  the  beam  should 
illuminate  a  large  amount  of  the  aperture.   Thus,  using  the  same 
system  as  one  would  use  when  taking  atmospheric  data,  a  recording 
of  the  system  alone  could  be  made.   This  would  be  most  helpful  in 
analyzing  single  photographs,  and  sequences  would  not  be  necessary. 

A  second  suggestion  is  to  convert  the  video  output  of  the 

television  camera  to  digital  data  for  later  machine  analysis.   The 

basic  equipment  is  available,  and  only  a  gating  circuit  and  sam- 

logic 
pling  sequence/need  be  added.   The  master  clock  that  generates 

the  pulse  train  used  to  synchronize  the  television  camera  runs  at 
110,25°  Hz;  the  first  frequency  used  in  the  television  synchroni- 
zation is  at  1575°  Hz,   The  sampling  pulses  could  be  selected  to 
run  off  this  master  clock  with  some  flexibility  as  to  the  desired 
frequency,  and  synchronization  with  the  television  camera  could 
still  be  retained.   With  the  data  in  digital  form  on  tape,  a  di- 
gital computer  could  be  used  to  analyze  the  data,  allowing  a  de- 
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tailed  analysis  over  the  entire  aperture. 

As  longer  paths  are  used,  the  problem  of  reflector  induced 
tilt  would  be  lessened.  The  increase  in  beam  diameter  would  not 
be  too  large  as  path  length  is  extended,  but  beam  motion  due  to 
reflector  tilts  would  be  more  apparent ,  The  motion  of  the  beam 
could  be  observed  and  data  could  be  taken  only  when  the  beam  is 
steady  on  the  target.   Thus  mirror  induced  tilt  would  be  minor. 
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APPENDIX  A 

This  appendix  contains  an  analysis  of  the  heterodyning  sys- 
tem in  more  detail  than  in  Chapter  2,  a  fuller  explanation  of  the 
data  analysis  system,  and  electronic  circuit  diagrams. 
A_^__l_Heterpdyne .  System  Analysis 

The  scheme  outlined  in  Figure  2  for  creating  a  spatial  in- 
terference pattern  proportional  to  the  phase  distortion  process 
caused  by  atmospheric  turbulence  can  be  analyzed  using  simple 
models. 

The  front  end  of  the  laser  (typically  the  output  of  the 
Spectra-Physics  Model  119  laser  would  be  0.1  milliwatt  at  the 
front  end  and  10  microwatts  at  the  back  end.)  is  collimated  and 
transmitted  through  the  turbulent  atmosphere  and  then  reflected 
back  over  the  same  path  by  the  down  range  reflector.   After  col- 
lection by  the  receiving  telescope  this  signal  is  mixed  with  a 
local  oscillator  signal  taken  from  the  back  end  of  the  laser. 
These  two  waves  can  be  approximately  described  before  transmission 
as  plane  waves  of  finite  spatial  extent.   As  the  finite  spatial 
extent  of  the  fields  should  not  cause  important  effects  for  this 
experiment,  it  is  not  explicitly  included.   Thus  in  complex  nota- 
tion 

I    ^Trans  =  %e      e  (A~1} 

-JB.7>t  (    } 

ILo     =  H2e        e 
where 

2tt   o) 


X  =  6328JI 
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The  atmospheric  effects  on  the  propagating  beam  include  dif- 
fraction, scattering,  and  absorption  losses;  atmospheric  turbulence 
modulation;  beam  spreading  due  to  imperfect  reflector-mirror  flat- 
ness; phase  modulations  due  to  mirror  motions  caused  by  vibra- 
tion.  Collection  of  the  returned  beam  will  introduce  further 
losses  due  to  the  finite  collecting  aperture  which  is  smaller  than 
the  beam  size.   Background  radiation  will  also  be  collected  in  the 
receiving  telescope.   The  various  effects  will  be  described  as  fol- 
lows '   _ 

Y(r,t) 

e        represents  the  turbulent  channel 

where 

Y(r,t)  =  X(r,t)  +  J<|>(r,t) 

X(r,t)   represents  the  amplitude  scintillation  process 
q)(r,t)   represents  the  phase  distortion  process 
L(r,t)   represents  losses  due  to 

1.  Diffraction 

2.  Particle  scattering 

3.  Molecular  absorption 

4.  Mirror  imperfections 

5.  Collection  losses 


e 


j*(r,t) 


iA  f§P^^-sen^s  phase  modulations  due  to  mirror  movements 
n 


E^CrjtJe        represents  the  wide  bandwidth  Gaussian  noise 

collected  by  the  receiver 


Islng  the  channel  model  of  Chapter  1,  the  transmitted  field  will  be 

received  as  _  _ 

jfyr   J  out   Y(r,t)    M(r,t) 
ER      (r,t)   =  L(r,t)lLe  e_     e  e 

ReC  J1        J€fc(r,t)    Jwt 

+  En(r,t)  e        e  U~3) 

The  collecting  telescope  will  minify  the  portion  of  this  re- 


oeived  beam  within  its  aperture.   As  long  as  one  looks  "close 
enough"  to  the  eyepiece,  the  minified,  image  is  a  clean  represen- 
tation of  illumination  in  the  object  space.   (Going  further  away 
gets  into  the  Fresnel  zone  of  the  secondary  mirror  and  aperture 
edge  effects.)   The  collected  field  is  mixed  with  a  local  oscil- 
lator field  by  a  beam  splitter  end  projected  on  a  plane  screen  for 
detection.   The  geometry  of  this  plane  screen  (Fig.  1)  is  selected 
so  that  z~0  on  the  surface;  the  transmitted  wave  approaches  the 
surface  normally  and  the  local  oscillator  signal  approaches  the 
surface  obliquely,  and  thus  has  a  direction  cosine  on  the  surface. 
The  fields  using  the  geometry  of  Figure  1  are 

ELo(p,t)  =  U2e     e  (A-k) 

-      _     y(pVO  jwt 

~Trans^p,t)  =  L(Pst)'aII1e      e 

J0n(p,t)  jiot 
+  |iEn(p,t)o   "     e  (A-5) 

Stotal  =  Elo^*)  +  £Trans(P,t)  (A-6) 

where 

p"  -   the  two-dimensional  vector  in  the  plane  of  the 
detector  surface 

\i  -   the  collecting  telescope-eyepiece  magnification 
The  detector,  in  this  case  film,  will  measure  the  intensity  which 

is  Hi, 

T   ,-  .  .    *Etotal^Etotal* 
Iin(p5t)  =  —————— 

Further,  a  black  box  model  of  this  detection  process  includes  a 
time  averaging  over  the  period  of  exposure  as  shown  in  Fig  A-l. 
T  is  the  period  of  the  shutter  opening  and  is  on  the  order  of  2«10 
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u 


EjP.t) 


n 


2L(p,t)u  U^Cp.tOe 


X(p,t) 


cos[©n(p,t) 


4»(p,t)-  ^(p,t)i 


and   2uU2En(p,t)cos[3p.  p  -©n(p,t)] 

The  integral  of  the  first  term,  the  square  of  the  envelope 
of  a  complex  Gaussian  process,  is  a  relatively  long-term  integra- 
tion of  a  very  wide  bandwidth  process;  thus  the  time  integration 
can  be  replaced  with  the  ensemble  average 

•T 


iJ/flcfctJlNKs     EJ>|e.<£0|1]  =    ff 


2 

where  K   is  proportional  to  the  average  power  in  the  relevant 

background  radiation. 

Similarly,  the  time  integrals  of  the  terms  containing  both 
E  (p,t)  and  ^„(p,t)  may  be  replaced  with  ensemble  averages  since 


n 


n 


both  E  (p,t)  and©jfp,t)  are  very  wide  bandwidth  processes  com- 


n 


pared    to   the   integration  time   T.      Thus,    for  example 

~r  f   ZJjl  Uz'Encfft)CDS[^f-o«Cf^]cJt    =      f  J  ^ik&ifa  &&[}#- 0*&$] 

Further,  the  two  processes  E  (p",t)  and©  (p,t)  are  independent, 


so  that 

E 


[*jL&E«<ffiQ*q,1-e«t$*  Ef[*A^E^]  E[^[J»f-G^4 


Also,(~)  (p",t)  is  reasonably  modeled  as  approximately  uniformly 
distributed  between  0  and  2rr  at  any  time  during  the  period  0  to  T. 
Thus  E  1  cos[B  •  p  -  Q   (p,t)]\  tends  towards  zero  and  the  time  inte- 
gral is  approximately  zero  over  the  entire  surface.   Similarly,  the 


>*cfr^ 


ZLCftt)^   U^Cf.-Oe  f'  W[©^(\*&--4><£t}-<K£t*] 


-57- 
terra  tends  towards  zero. 

This  integration  model  of  the  detector  implies  that  the  de- 
tector records  energy.   For  this  model,  since  the  processes  are 
slow  compared  to  the  shutter  time,  a  typical  intensity,  I(p~),  has 
been  used  with  the  understanding  that  multiplication  by  T,  the 
integration  time  constant,  will  yield  the  proper  energy  terms. 

The  resulting  output  from  this  integration  model  is  thus 

if,,2     22_2  2X(p~)    2 
I(P)  =   27)     \     \H2\      +   M  L  (p)l~l\  e      +  N 

+   2ML(p)    luJlUgje  cos[Bp.p   -  (j?{  p)-  ^(  p)    ]V      (A-ll) 

The  loss  term  L(p)  is  reasonably  modeled  as  a  uniform  quan- 
tity over  the  strong  portion  of  the  transmitted  beam  since  scat- 
tering and  molecular  absorption  are  essentially  non-preferential 
and  are  less  apparent  in  the  strong  portion  of  the  beam.   As  it 
is   this  portion  of  the  beam  that  will  be  collected,  the  term  L( p ) 
will  be  replaced  with  L. 

The  term  ty(p)  represents  the  phase  modulation  due  to  mirror 
vibrations  further  acted  on  by  the  turbulent  channel.   Due  to 
mirror  mount  massiveness,  the  value  of  ty(p»t)  is  constant  with 
time  over  the  shuttering  period.   The  spatial  variation  of  ty(p) 
is  critical,  however,  since  the  mirror  movements  that  can  occur 
cause  a  tilt  on  the  reflected  wave.   This  tilt  is  of  the  same  form 
that  one  expects  as  the  first  order  effect  of  atmosphere  dis- 
tortion. 

There  are     arguments  that  seem  to  indicate  that  while  this 
tilt  caused  by  the  mirror  is  a  limitation  on  the  system,  the  mea- 


purements  of  the  atmospheric  effects  can  still  be  made.   The  first 
is  that  typical  tilt  changes  between  photographs  are  on  the  order 
of  15  microradians.   If  this  tilt  change  were  due  only  to  mirror 
movements,  this  would  mean  a  beam  movement  of  0.  05  inch  at  the 
telescope.   The  overlap  of  the  returned,  beam  over  the  receiving 
was  approximately  0. 5  inches,  meaning  that  this  movement  could, 
not  be  seen.   However,  more  significant  tilt  changes  would  be  100 
plus  microradians,  which  translates  to  a  0. 5  inch  travel  at  the 
telescope.   Thus  this  movement  would  be  seen  and  neglected.   Thus 
'it  is  recognized  that  mirror  motions  are  critical,  but  for  this 
path  do  not  swamp  out  larger  values  of  atmospheric  tilt.   There- 
fore, ty(p)  will  be  lumped  in(j>(p"),  and  the  intensity  can  be  written 

as 

2  2  2    2,2    2X(p) 


2?|  " 


L  Ull    e 

X(p) 


+   2uLluJluJe  cos[8p.  p   -^>(p)-a]|  (A-12) 

where  the  intensity  contains  a  term  phase  modulated  by  the  atmos- 
pheric phase  distortion. 
A  -  2  Electronics 

The  data  analysis  system  was  built  around  the  Packard  Bell 
920  television  camera  and  is  shown  in  block  diagram  form  in  Fig. 
A-2.   As  produced.,  the  camera  has  a  211   random  interlace  which 
was  altered  to  a  single  vertical  scan  to  allow  accurate  reading 
of  the  fringe  pattern.   The  sweep  of  the  television  camera  was 
made  into  256  horizontal  lines  each  vertical  sweep,  with  synchron- 
ization between  the  horizontal  and  vertical  rasters.   The  circuits 
used  for  this  system  were  built  around  RTL  digital  integrated- 
clrcuits. 
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Figure  A-?.  Data  Analysis  System  Bloclc  Diagram 
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A   master  cloak  shown  in  Fig,  A-3  was  used  to  generate  the 
pulse  train  for  the  system. 

A   series  of  counters  and  pulse  shapers  were  used  to  generate 
synchronization  pulse  trains  from  the  master  clock  signal.   In 
each  case  a  synchronous  counter  was  used  to  minimize  countdown 
delays  inherent  in  sequential  counters.   Figs.  A-'+,  i\-5  and  A-6 
are  circuit  diagrams  for  these  counters  and  pulse  shapers.   In 
each  of  the  circuits  the  counting  is  done  by  the  JK  flip-flops 
gated  on  and  off  by  the  logic.   The  shapers  were  pulse  generators 
with  a  variable  resistor  used  to  set  the  output  pulse  width. 

The  synchronization  pulses  were  put  into  Schmidt  trigger 
circuits  to  provide  some  power  gain  and  also  to  buffer  the  count- 
er circuits. 

Using  the  system  outlined  above,  the  television  camera  scan 
<  was  made  to  cover  the  256  horizontal  lines  in  each  vertical  sweep 
I with  no  interlace.   The  starting  point  for  each  vertical  sweep 
was  the  same  so  that  the  video  output  could  be  displayed  on  an 
oscilloscope,  using  the  synchronization  pulses  as  triggers,  with 
precision  and  stability. 

Before  display,  the  video  output  of  the  television  camera 
was  put  into  a  low  pass  filter  with  a  3  db  point  at  k   MHz.   The 
low  pass  filtering  was  done  to  remove  a  large  amount  of  high 
frequency  hash.   Fig.  A-?  shows  the  improvement  of  the  display  pro- 
duced by  this  filtering. 


OFFTiMC   £   ^f^ 

Figure  A- 3  Master  Clock  and.  Pulse  Generator  for 
Synchronization  System 
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Figure  A-iJ.     Eight   Counter  and   Pulse  Generator 
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Figure  A-&    Seven  Counter  and  Pulse  Generator 


OOTPcJT 


Figure  A~^3  Four  Counter  and  Pulse  Generator 
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Figure  A- 7. 


Improvement  in  Picture  Quality  from  Insertion 
of  Low  Pass  Filter  at  the  Television  Video 
Output.   Upper  Picture  is  without  filtering; 
Lower  Picture  is  with  filtering. 
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This  appendix  contains  a  discussion  of  the  data  analysis  sys- 
tem limitations,, 
B  -  1  Fringe  Pattern  Model 


As  shown  in  Appendix  A,  the  intensity  for  plane  wave  models 
is 


,  r  2  2    2  2    2  2x(p) 

I(P)  -^{llfel   +"»   +^}1%I 


+  2WLJUJ  JU2|e    cos[Bp-  p  -<j>(p)-a]j     (A-12) 
The  fringes  result  from  the  spatial  sinusoid,  which  is  proportional 
to  the  product  of  the  two  signals.   The  fringes  ride  a  bias  which 
is  the  sum  of  the  intensities  of  the  two  signals.   The  terms  re- 
presenting atmospheric  turbulence  effects  can  be  thought  of  as 
minified  representations  of  the  processes  at  the  objective  of  the 
receiving  telescope.   Thus 
X(p)  =  X(p'/u) 
(j>(p)  ^(P'AO 
where 

p   is  the  two  dimensional  vector  at  the  detector  surface 

p'   is  the  two  dimensional  vector  at  the  objective  plane 
oriented  so  that  the  geometry,  when  propagated  through 
the  receiver  optics,  matches  that  of  p. 

|j  represents  the  magnification  of  the  telescope-eyepiece 
combination 

Further,  aligning  the  geometry  of  this  model  so  that  the  fringe 

lines  coincide  with  one  of  the  axes  in  the  p  plane  simplifies  the 

notation.   With  this  assumption,  the  intensity  is 
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r        2    2 

1' 


,  r   ,2     2     2  2     2  2X(p) 
I(p>  =  2^{|-Ull   +  N   +  »  L  \ty    6 


where 


+  2uL^U1J|U2lcos[8xx  -(j)(p)~a] 

Bx  =  lB|(Bx/lBP  =  |Blsin£=  jpsinl*  2nlA 

£  represents  the  angle  between  the  two  interfering  waves 
and  is  set  by  the  experimenter  so  that  £<£<! 

By  looking  at  a  very  narrow  strip  in  the  y  dimension  so  that  vari- 
ations in  y  do  not  have  a  chance  to  occur,  the  intensity  model  can 
be  written  as 

_  /-        2  2  2   2  2   2X(xyy    ) 

I(P)   =  2»j|Hi|      +  N     +  M  L  |U2|    e 

.     ,  X(x,y   )  p    0  ? 

+  2ULJUJ  [ujje  °  cos[-2^-  -  cj>(x,yo)-a]  j         (B-l) 

The  problem  of  seeing  the  fringe  pattern  is  related  to  the 

depth  of  modulation,  which  can  be  shown  in  this  model  to  depend  on 

the  visibility  ratio  V,  defined  as 

V  _  amplitude  of  sinusoid 

amplitude  of  bias 

UJ  +  N  +  u  £  I  Uj  e 
For  this  experiment,  N  was  negligible  compared  to  the  other  terms 

in  the  denominator  and  will  be  neglected.   The  other  terms  in  the 
denominator  represent  the  local  oscillator  and  received  beam  in- 
tensities.  The  shadow  pattern  of  these  intensities  can  be  seen 
during  experimentation.   From  observation,  the  shadow  pattern  of 
the  received  field  is  brighter  than  the  local  oscillator  by  an 
estimated  factor  of  three  to  four.   On  a  longer  path,  losses  will 
surely  be  more  severe  and  this  ratio  would  be  reversed.   If  Eq„  B~2 
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is  manipulated  slightly,  and  the  information  on  the  ratio  of 
the  two  beam  intensities  incorporated,  the  visibility  ratio  be- 
comes 

X(p) 


V  =   2    (uLJuJ  /|U2|  )e 


2.2 


1   -I-   M   L^juJ  C      2X(p) 

.X<P>  =        4eX(p) 


2  J£e/ 

—  Kr2XfpT  1  +4e2X(p) 


(B-3) 


It  can  be  seen  that  the  scintillation  process  is  the  fundamental 
limitation  on  modulation  depth  for  this  path.   This  is  plotted  in 
Fig.  B-l. 

V 
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Figure  B-l   Plot  of  Visibility  Ratio  with  changes  in 
the  value  of  scintillation 
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A   further  limitation  peculiar  to  this  particular  data 
analysis  system  was  the  resolution  available  using  an  oscillo- 
scope display  for  detailed  analysis.   The  scope  face  was  ten 
centimeters  across  and  distinctive  portions  of  the  waveform 
could  be  located  accurately  to  one-half  millimeter  as  long  as 
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the  surrounding  clutter  was  not  too  bad.   However,  at  least  a 

three  to  five  millimeter  separation  was  required  before  it  was 
possible  to  reliably  distinguish  peaks  from  clutter  for  the  typi- 
cal depths  of  modulation  produced  by  the  scanned  fringe  lines. 
This  means  that  20  to  33  lines  maximum  could  be  handled  on  the 
scope.   It  was  found  that  (see  section  3~1 )  the  effects  of  the 
atmospheric  turbulence  on  the  fringe  perturbations  did  not  become 
easily  visible  until  the  number  of  fringes  across  the  aperture 
became  twenty  or  less.   Thus  the  use  of  the  scope  was  not  a 
serious  limitation  for  this  experiment. 
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